NASA 


CO 


TECHNICAL 

REPORT 



m NASA TR R-310 


tO'AN COPV: I . 
AFWL <V Qi 
KiRTLAND A!°j 

03 i 
Ln 

a 

nj 


Jb 


2 * 
» I 




n 

CD 

z 

2 


A SEMIEMPIRICAL COLLISION 
MODEL FOR PLASMAS 


by Willard E. Meador 

Langley Research Center 
Langley Station , Hampton , Va. 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION % WASHINGTON, D. C. * MAY 1969 



TECH LIBRARY KAFB, NM 



DDbflSDa 


>■« v J. V 


A SEMIEMPIRIC AL COLLISION MODEL FOR PLASMAS 

By Willard E. Meador 

Langley Research Center 
Langley Station, Hampton, Va. 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


For sale by the Clearinghouse for Federal Scientific and Technical Information 
Springfield, Virginia 22151 — CFSTI price $3.00 



CONTENTS 


SUMMARY 1 

INTRODUCTION 1 

SYMBOLS 3 

THE COLLISION MODEL 7 

Derivation of Model 7 

Determination of K and £j 9 

A NEW COLLISION TIME 10 

THE MAGNETIC PROBLEM 14 

TEMPERATURE EFFECTS 17 

Thermal Conduction 17 

Thermal Diffusion 18 

ALTERNATE DETERMINATION OF £ 19 

CONCLUDING REMARKS 20 

APPENDIX A - COLLISION INTEGRALS 21 

APPENDIX B - THE KROOK COLLISION MODEL 23 

APPENDIX C - CHAPMAN-ENSKOG SOLUTIONS 26 

APPENDIX D - THE ONSAGER MAXIMIZATION PRINCIPLE 29 

APPENDIX E - THE MEAN-FREE-PATH MODEL 30 

REFERENCES 32 


iii 



A SEMIEMPIRICAL COLLISION MODEL FOR PLASMAS 


By Willard E. Meador 
Langley Research Center 

SUMMARY 

A collision model similar in many respects to that of a Lorentz gas is developed 
for plasmas. Two assumptions are made: (1) the totality of interparticle force laws 

(including electron-electron) can be replaced with a single parameter determined by the 
unique ratio of generalized electric field to temperature gradient that can exist in a 
bounded, nonmagnetic, and currentless plasma of given composition; (2) the electron 
dynamics governed by such an average potential are Lorentz-like in the sense that the 
effective electron collision partners behave as heavy scattering centers. With the col- 
lision model thus specified and held fixed, the first-order kinetic equation is solved 
exactly for a succession of more complex plasma conditions; the corresponding distribu- 
tion functions are used to predict such diverse properties as entropy, entropy and diffu- 
sion relaxation times, electromagnetic tensor conductivities, thermal conductivity, and 
thermal diffusion. Since the results are generally comparable to the third Sonine approx- 
imation and nearly always better than the second approximation (thermal conductivity 
excepted) for a fully ionized gas, the method apparently provides an adequate scheme for 
utilizing experimental data when microscopic collision parameters are scarcely known 
(e.g., electron — neutral-atom forces) and for directly extending previous calculations to 
more complex conditions. Interesting sidelights include the discovery that the elementary 
mean-free-path formulation of the Hall conductivity is left intact for small magnetic fields 
if the familiar collision time is replaced by a new definition associated with entropy 
relaxation. 


INTRODUCTION 

In the past few years, many authors have discussed statistical models as approxi- 
mations to Boltzmann collision integrals for solving problems in rarefied gas dynamics 
(e.g., ref. 1). The philosophy is summarized briefly as follows: Whereas the exact col- 

lision integrals treat in detail the geometry of every binary encounter, the essential fea- 
tures can be described almost as well by first taking a statistical average over these 
encounters. Aside from providing validity tests for various approximate distribution 
functions, the models could prove useful in their own right for such purposes as adapting 
kinetic theory to experimental data in lieu of unknown microscopic parameters, directly 



extending previous calculations to more complex gas conditions (e.g., the extension of 
nonmagnetic results to electromagnetic fields), and solving problems that do not respond 
to classical perturbation methods. In addition, the flux expressions derived from simple 
models are often far more concise than those from polynomial expansions and are there- 
fore easier to interpret. 

Although some of the concepts can be generalized to neutral gases, the present 
investigation is concerned entirely with plasmas and particularly with those processes 
controlled by electrons. The development of a simple collision model and therefrom a 
reliable semiempirical electron distribution function is motivated by the following con- 
siderations: 

(1) Electron-electron encounters prevent closed-form solutions of the Boltzmann 
equation and thus obscure important characteristics that might be better understood in a 
more empirical scheme. 

(2) Most of the electron — neutral-atom interaction potentials are scarcely known. 

However, one must recognize that the experimental data readily accessible from plasmas 
are very limited and, consequently, questions arise as to the type of data, the amount of 
data necessary for given purposes, and the theoretical framework in which the data are 
best represented. Information theory (ref. 1), for example, maximizes the observer’s 
uncertainty about all plasma properties except those actually measured, but the entropy 
and distribution function, and so forth, so determined bear little resemblance to the true 
quantities if the experimental constraints are inadequate. Besides the usual conservation 
requirements, prior knowledge is ultimately assumed in reference 1 of all second-order 
velocity moments (including off-diagonal elements) in a multicomponent mixture and also 
the interparticle force laws for the specification, though somewhat arbitrarily, of certain 
other parameters. Such methods are of little value when a purpose of the semiempirical 
procedure is to avoid unknown potentials. 

The problem then is to enhance the adequacy of a small number of practical con- 
straints, at least one of which supplants theoretical cross sections, by somehow intro- 
ducing a bias in the observer's attitude toward things not measured. Some indication of 
the proper bias can be deduced from the Boltzmann equation if a sufficiently simple col- 
lision model is adopted. The intent of the present research is to develop such a model 
based on the following assumptions: (1) the totality of interparticle force laws (including 

electron-electron) can be replaced with a single parameter determined by the unique ratio 
of generalized electric field to temperature gradient that can exist in a bounded, nonmag- 
netic, and currentless plasma of given composition; (2) the electron dynamics governed by 
this average potential are Lorentz-like in the sense that the effective electron collision 
partners behave as heavy scattering centers. With the collision model thus specified and 
held fixed, the first-order Chapman-Enskog equation is solved exactly for a succession 
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of more complex plasma conditions; the corresponding distribution functions are used to 
predict such diverse quantities as entropy, entropy and diffusion relaxation times, elec- 
tromagnetic tensor conductivities, thermal conductivity, and thermal diffusion. 

Since the results are generally comparable to the third Sonine approximation and 
nearly always better than the second approximation (thermal conductivity excepted) for a 
fully ionized gas, the method apparently provides an adequate scheme for utilizing exper- 
imental data in multicomponent partially ionized plasmas and for directly extending pre- 
vious calculations to more complex conditions. Interesting sidelights include the dis- 
covery that the elementary mean-free-path formulation of the Hall conductivity is left 
intact for small magnetic fields if the familiar collision time is replaced by a new defi- 
nition associated with entropy relaxation. 

For simplicity, it is assumed herein that heavy particles (e.g. , ions and neutral 
atoms) are infinitely massive and are at rest relative to the laboratory. 

SYMBOLS 

a parameter such that a = 0 refers to a Lorentz fully ionized plasma and 

a = 2 1/2 refers to a real fully ionized plasma 

b impact parameter 

B magnetic field 

B unit vector in direction of magnetic field 

c particle velocity 

Di,D2 ? Dg functions defined in equation (C2) 
e negative of electron charge 

E electric field 

f distribution function 

f(°) Maxwellian distribution function 

g* function defined in equation (5) 
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enthalpy flux 


unit vector in x-, y-, and z-direction, respectively, of Cartesian coordinate 
system 

integral defined by equation (C8) 
electron current density 
Boltzmann's constant 

parameter in the collision model defined by equation (Cl) 

particle mass 

particle number density 

electron partial pressure 

kinetic energy flux 

-2 

combination Rjg 1*04-^22 ^ij integrals 

limit of Ijj integral as co approaches zero 

entropy density 

equilibrium entropy density 

collisional production rate of entropy density 

time 

electron temperature 

measure of average energy with which electrons emerge from encounters 
with effective collision partners j 


u e j 


x 


0 

/3ej 


r 


y E 


7t 

6 

6 E 


6 


T 


€ 


e 


€ H 

T 1 

e d 


e t 


e 


average velocity with which electrons emerge from encounters with effective 
collision partners j 

electron diffusion velocity 

integration variable 

reduced electron diffusion velocity 

reduced average velocity with which electrons emerge from encounters with 
effective collision partners j 

reduced electron particle velocity 

Spitzer-Harm reduction factor for scalar electrical conductivity 
Spitzer-Harm reduction factor for thermal diffusion 
variational operator 

Spitzer-Harm reduction factor for kinetic energy flux due to electric field 
and/or pressure gradient 

Spitzer-Harm reduction factor for thermal conductivity 
azimuthal angle in collision dynamics 

effective electron driving force (generalized electric field) defined by 
equation (C4) 

component of ef parallel to magnetic field 
component of T perpendicular to magnetic field 
electron diffusion energy density 
total nonrandom electron energy density 
polar angle in spherical coordinates 
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4 

a 

a H 

a l 

°H 

a L (~) 


<P 

<Pe 

X 

co 

Subscripts: 

c 

e 


Lagrangian multiplier 

effective interaction parameter for electron collisions 

scalar electrical conductivity 

parallel (to magnetic field) conductivity, <7|| = o 

perpendicular (to magnetic field) conductivity 

Hall conductivity 

exact first-order conductivity (infinite Sonine approximation) for a Lorentz 
gas 

collision time equal to t ct ( 1), that is, first Sonine approximation to t g 

collision time based on entropy relaxation 

collision time deduced from scalar electrical conductivity 

azimuthal angle in spherical coordinates 

electron first-order perturbation function 

scattering angle 

cyclotron frequency for electrons 

specifies effect due to collisions 

electrons 

index 


6 



j 


index; also particles of type j 


L Lorentz gas, that is, no electron-electron collisions 

Arguments 1, 2, 3, . . ., °° refer to level of Sonine approximation. 

Primes with symbols refer to quantities after a collision. 

The notation <A> refers to average of quantity A over velocity space. 

THE COLLISION MODEL 

Derivation of Model 

The Boltzmann collision integral (ref. 2) for a two-component Lorentz plasma (no 
electron-electron encounters) is written to first order and for fixed scattering centers j 
as 


(inr) ej = J ( f eV - f e f j)| ? e ~ Sj|b db de dcj 
= ^(fe* - fe)fjc e b db de dcj 

= n )^ <0) (^ £ ) 1/2 I(»’e , -^)»dbd e 


(1) 


where primed symbols refer to quantities after a collision (as opposed to unprimed sym- 
bols for quantities before a collision), f e ^ is the Maxwellian distribution function for 
electrons given by 


(o) _ 


m. 


,3/2 


= n f 


-r 2 


( 27 rkT e . 


( 2 ) 


Y is the reduced electron particle velocity expressed as 


y = 



and cp e is the electron perturbation function defined by 


(3) 


(4) 


f e = f e^ (l + <Pe) 
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Equation (1) is further simplified by observing that cp e is always dependent upon 
y in the following manner (ref. 3) for the assumptions stated in the last paragraph of the 
Introduction: 

cp e = g(y , fields and gradients) • y (5) 

Accordingly, 



if | is the power in the electron — heavy-particle inverse power interaction potential. 
(See appendix A.) The parameter K embraces a collection of thermodynamic param- 
eters and also may be a function of ij. 

No matter how many heavy components the plasma contains, the assumption that £ 
can be so adjusted as to make the far right side of equation (6) a fair representation of 
the sum of all electron — heavy-particle collision integrals appears justified. It seems 
less likely, on the other hand, because of the more complex integrals, that electron- 
electron encounters can be thus included. Nevertheless, 



= -K y 


l-4/« f (o) 


<Pt 


(7) 


is adopted as the complete (j = e included) real-gas collision model, the validity of which 
will be tested by macroscopic comparisons. This collision model definitely establishes 
K and £ as empirical parameters and has an advantage over other simple models in 
that £ is immediately identifiable and (di e jdtj^ is the exact first-order expression 

if the plasma is truly Lorentz with particles of only a single heavy species serving as 
active electron scattering centers. This condition holds true either for a slightly ionized 
pure gas (negligible electron-electron and electron-ion effects) or for full ionization with 
large ionic charge. 

It is further evident that equation (7) reduces to the well-known Krook model (cf., 
ref. 4 and appendix B) only when £ = 4. The present method is thus a generalization of 
Krook's work and shows his model to be restricted to Maxwellian molecules (£ = 4) for a 
Lorentz plasma and to the first Sonine approximation to the Chapman-Enskog solution 
(ref. 2 and, also, appendix B) for a real gas. Since the first Sonine approximation is 
completely inadequate in its description of many plasma features, the additional y- 
dependent factor in equation (7) should constitute a major improvement. 
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Determination of K and £ 


Appendix C provides the exact solution of the Chapman-Enskog first-order perturba- 
tion equation corresponding to the collision model of equation (7) and to the application of 
a pressure gradient and an electromagnetic field to a zero temperature gradient plasma. 

Of more immediate interest, however, is the nonmagnetic plasma, the perturbation func- 
tion for which is found by setting oo equal to zero in equation (C9) and using Ohm's law 
j = -en e Vg = oT to write 


<Pe = 2R 04 R 13 _ • y 


( 8 ) 


The Rjj integrals are defined in appendix C and 


/3 = 



(9) 


is the reduced electron diffusion velocity. 

In addition, K is determined for equation (C7) by the requirement that the velocity 
moment of cp e , which originally depended on K, must yield the correct electron diffusion. 
Therefore, from equation (7) 



R 13 

T a^04 


y 


l-4/if (O) 


<Pt 


(10) 


where r g is the familiar collision time defined in terms of the electrical conductivity 


a = e^n e T a ym e (11) 

Only £ remains unspecified and, therefore, the entire burden of accounting for the aver- 
age interparticle collision properties falls on that parameter. 

It seems reasonable to propose that an adequate value of £ is obtained by requiring 
the energy moment of equation (8) to duplicate the experimental kinetic energy flux q. 

One therefore computes 

q cc (^2 - n e -1 Jy 2 >5 e ^ 0 Ve d ^e = 2 £ -1 (l + t)P (12) 

the energy-flux — enthalpy-flux difference 

(y 2 r) - 1 <f> = (y 2 r) - f fi = (4 (is) 
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for future reference, and finally 

2r 1 £= (y 2 r) - 2/s d 4 ) 

from either equation (12) or equation (13). Each expression refers explicitly to the non- 
magnetic, zero-temperature-gradient gas. 

Equations (10) and (14) with n e , T e , and measurements of 0 and repre- 

sent the complete semiempirical collision model. Although equation (14) obviously does 
not refer to the bounded, nonmagnetic, and currentless plasma postulated in the Introduc- 
tion for the determination of |, the two sets of conditions are subsequently shown to be 
exactly equivalent in that respect — that is, the values found for £ are identical in the 
two situations. Also, £ may not be the power in any particular inverse power interac- 
tion potential, nor is it necessary to suppose that the constituent force laws are of that 
type. Macroscopic data have replaced the microscopic dynamics. 

As mentioned previously, the special plasma conditions underlying equation (14) 
serve only to fix experimentally the numerical value of ij, which is constant for a given 
set of plasma constituents. The electron first-order perturbation function cp e therefore 
remains the lone condition-dependent function in equation (10), and it must be determined 
separately for each problem considered. The procedure is always the same: 

Replace the Boltzmann collision integrals with equation (10), substitute the appropri- 
ate fields and gradients into the first-order kinetic equation, and solve the resulting 
expression exactly. 

Illustrations and tests of this theory are provided herein. 

A NEW COLLISION TIME 


As a preliminary to the first application and test of the above semiempirical theory, 
it is instructive to compute the rate at which collisions produce entropy density s c ; 
accordingly, from equations (8) and (10) 


*. ■ - J-0 * 4i), »» 

The rate of dissipation of the diffusion energy density e d is then related to that energy by 


\ T a T e®c = e d = 1 n e m e v e 2 


(16) 


It thus appears that r a measures quite adequately the relaxation time scale for 
entropy, at least to the extent that e d represents the complete nonrandom electron 
energy density. But it will now be shown that this condition does not hold. The best way 
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of defining the total nonrandom electron energy density e^. is through the following rela- 
tion involving the difference between the entropy density s and its equilibrium value 
s(°): 

f e ln(f e 


e t = T e (s (o > - s) = kT e J 




f e ( °V e 2 d? e = ^d 


(17) 


The far right side of equation (17) is obtained from the substitution of equation (8) for cp 

-2 e 

and R for Rj 3 ~ ^04^22* (See appendix C.) 

A convenient connection between equations (16) and (17) and the concepts of irrevers- 
ible thermodynamics is provided by Onsager T s fundamental principle, which states that the 
collisional production rate of entropy density is a maximum (ref. 5, pp. 116 and 145-148). 
See also appendix D where it is shown that the present distribution function can be obtained 
by the maximization of s c subject to the constraint that /3 is known. This principle 
suggests a new collision time in equation (16) of the form 

T s = e d" le t T a = RT a ( 18 ) 

so that the entire is subject to dissipation in the manner 

\ T s T e s c = e t (19) 


Since T s/ T a excee ds unity (fig. 1) for all but Maxwellian molecules (£ =4), 
entropy equilibrization generally must lag that of electron diffusion. Such behavior is not 
surprising in view of the many collisions normally required to produce randomization. 

E. A. Mason of Brown University has suggested the possibility of r g - t q being a cor- 
rection for the so-called persistence of velocities in free-flight or free-path theory. If 
this suggestion is valid, the nonrandom energy might be explained in terms of the 

,f memory M of an electron, that is, the time required for an electron to forget its original 
speed and direction. In any event, the method for calculating r s is well prescribed by 
the relation 


and the parameter is 
subsequently. 


TS 2T e( s(0) - S J~ kT ^j* fe (o)^ e 2 d - 


n e m e v e 2 n e m e v e 2 


shown to occur naturally in the Hall conductivity discussed 


( 20 ) 
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Figure 1.- Ratio of collision times as a function of effective collision parameter. 


Meanwhile, the ratio t s j7 a can be computed and compared by several methods as 
a first test of the semiempirical distribution function. It is convenient for this purpose 
to express equation (14) in terms of the Spitzer-Harm relations (ref. 6) 


1 3 = 


m f 


en e \2kT e 


1/2 


J = - 


m c 


en P \2kT 


1/2 


( m 


ae = - 


en e \2kT e 


1/2 


a L ?= I'e^L 


( 21 ) 



( 22 ) 


between real and Lorentz plasmas operating under identical electric field and thermody- 
namic conditions. Hence, 
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2 r 1 r E ^ L = 6 E <r 2 f> L - 2 r E 0 L 


(23) 



Although the data necessary for experimental verification of equations (10) and (23) 
are not available, one can still check their accuracy on the basis of the Spitzer-Harm 
exact first-order solution for a fully ionized gas. The parameters y^, = 0.5816 and 
6g = 0.4652 are known for this plasma, as is the Lorentz quantity 

(r 2 r) L = 4j L (24) 


from equation (12) with | = 1; consequently, £ has the value 


- (We'^e “ 


1.6674 


(25) 


Electron-electron effects are therefore quite substantial in their adjustment of the origi- 
nal (£ = 1) electron-ion force law. This result implies that the subsequent tests of equa- 
tions (10) and (25) are indeed meaningful, even though the plasma is without neutral 
components. 

The integrals comprising R in equation (18) were computed numerically for the 
value of £ of equation (25) to yield Rq^ = 0.66467, R 13 = 1.2113, and R22 = 2.6461; 
accordingly, 

t ct -1t s = 1.199 (26) 


In like manner, substitution into equation (20) of equation (Cll) of appendix C with co = 0 
gives the second Sonine approximation 


T s(2) = 16a 2 + 104a + 259 _ j 25Q 
T a( 2 ) 16a 2 + 104a + 169 


(27) 


when a equals 2 4 / 2 which corresponds to the real plasma. Although the agreement is 
already reasonable, it is much improved by a similar comparison with the third Sonine 
approximation 


T s( 3 ) _ 20736a 4 + 173952a 3 + 518392a 2 + 370856a + 93169 
T CT (3) 20736a 4 + 173952a 3 + 427312a 2 + 262136a + 47089 


(28) 


obtained from the matrix elements of Schweitzer and Mitchner (ref. 3). 

The advantage of the present technique is even more pronounced for the Lorentz gas 
(i; = 1) for which equation (18) yields the exact first-order result 


13 



(29) 


T -1 T - 315ff 
s 512 


1.933 


as compared with the values 


T s( 2 ) 

T 0r( 2 ) 


1.533 


T s( 3 ) 

T CT (3) 


1.979 


(30) 


which correspond to a = 0 in equations (27) and (28). Equations (26) to (30) also illus- 
trate the large sensitivity of r s /T a to the effective interaction potential (ij = 1.6674 
and 1), as do the following ratios from equation (18) representing Maxwellian (£ =4), 
Coulombic (i; = 1), and rigid-sphere (£ = <*>) forces, respectively: 

T s/ T cr ~ I* 000 7 s/ 7 a = 1-933 t S/ /t ct = 1.178 (31) 

Similar conclusions apply as well to other plasma properties (e.g. , thermal and magnetic), 
so that the previously stated purpose of using experimental data to avoid scarcely known 
collision parameters indeed has merit. 

It is further observed that the first Sonine approximation always gives unity for 
T s/ T a> regardless of £, and so do equations (27) and (28) when a is very large. Since 
large values of a absolve the ions of their scattering responsibilities, such results sup- 
port the well-known claim that the first Sonine approximation is adequate for pure gases 
and for mixtures with small mass disparities. 

THE MAGNETIC PROBLEM 

A more exacting test of equations (10) and (14) or (25) is obtained by the application 
of an electromagnetic field. As explained earlier, the method will be useful for directly 
extending nonmagnetic measurements and previous nonmagnetic calculations to the full 
field if the specification of £ in terms of the nonmagnetic energy-flux — diffusion-velocity 
ratio is adequate. Again the fully ionized gas is considered, although one should recognize 
that other types of plasmas merely involve different empirical or theoretical values of yg 
and 5g in equation (25). 

The complete first-order solution of the electronic Lorentz-like Boltzmann equation 
for a plasma in an electromagnetic field is derived in appendix C, as is the current density 

T = oj|?j| + CTj?^ + ct h B xe = + r 13 _ 1 i 13 ~| + WT cr R 04 R 13 ~ 2l 22 B x?) (32) 
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Additional notation includes the parallel (to 6) , perpendicular, and Hall conductivities 
(cr||, a x , and Ojj, respectively^, the cyclotron frequency u> = eB^m e , and the integral 

j.. . p(41/SH(l ♦ »2 T( 2 Ro4 2 Rl3 -2 x <8A)-2) e - 2 (33) 


The limit of this integral, as cu approaches zero, is Rjj. One can also write the fol- 
lowing combination of equations (18) and (32): 

T = + R 13~ 1]t 13^L + wT s r 22 _1i 22B xTj (34) 

Since I^j and R^j are identical when quadratic and higher powers of cot 0 are 
neglected, equation (34) becomes 


j = cr(V 


+ 0JT S B X € 


(35) 


for small magnetic fields; likewise, the second So nine approximation of equations (27) 
and (C12) yields the same formal result. These results are interesting for at least two 
reasons: (1) the appearance of r s rather than r a in the Hall conductivity indicates a 

preference for the relaxation time of entropy over that of diffusion; (2) except for the new 
collision time, equation (35) is exactly Spitzer T s formula for the simple mean-free-path 
current density (ref. 7; see also appendix E). In addition, a comparison of the numerical 
results of equations (26) to (30) suggests that the present technique for computing to 

linear ojt s might be better than the third Sonine approximation. 

Although equation (34) still provides the most concise description, the quadratic and 
higher effects of large magnetic fields are difficult to interpret. Even so, the conductiv- 
ity ratios are found quite easily from equation (32) by using the nonmagnetic value of £ 
from equation (25) and y-g(°°) = 0.5816 in the following formulas: 


and 


g i 

a L (°°) 


- r E (°°)Ri3(4) 


(36) 


= r E ( oo )wT o R 0 4 R 13 (^)" 2 l22(^ WT a) (37) 

o r L(°°) 

where 0 ^( 0 °) is the exact Lorentz conductivity y^^cr = The results obtained 

with equations (36) and (37) are listed in the last column in tables I and n, respectively, 
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for the extension of Spitzer and Harm’s scalar conductivity (ref. 6) to the electromagnetic 
tensor conductivity at several values of a n. 

Also listed in tables I and II are calculations from the first Sonine expressions (see 
appendix E) 

^d) r E U) 

1 + 0)2 r a ( l)^ 
ct h ( 1) r E (l)o)r a (i) 

<r L < “> 1 + cA^l) 2 

and the second, third, and sixth approximations of Kaneko (ref. 8). Besides 7^(1) = 3 tt/ 32 
from reference 6, a few conversion factors are necessary: the r a (l) of equation (38) is 
equal to the r of tables I and II and each entry in table I of reference 8 should be multi- 
plied by ±yg(l), where the minus sign is used in the Hall conductivity. 

The comparisons are generally favorable and show equations (36) and (37) to be 
nearly always superior to the second Sonine approximation. One can thus proceed with 
some assurance to the extensions and semiempirical applications originally intended for 
the present method. 


(38a) 


(38b) 


TABLE L- PERPENDICULAR CONDUCTIVITIES 


TABLE II.- HALL CONDUCTIVITIES 


ojrt 

^(D/ulH 


<q(3)/a L^°°) 


CT i/°L ( " ) 
(eq. (36)) 

0.0 

0.295 

0.569 

0.574 

0.581 

0.582 

.2 

.283 

.458 

.468 

.471 

.471 

.5 

.236 

.236 

.253 

.256 

.266 

1.0 

.147 

.100 

.106 

.110 

.118 

2.0 

.0589 

.0406 

.0367 

.0381 

.0393 

4.0 ; 

.0173 

.0147 

.0126 

.0117 

.0111 

6.0 

.0080 

.0073 

.0066 

.0057 

.0051 


tThe first-order collision time r is related to r a by 
the ratio T 0 jr = 32(3t r)" Vgj where y £ is given in reference 6; 
in particular, y E (°°) yields Tg^r - 1.975 for the conversion of 
the independent variable in equation (36). 


cjrt 

<T H <y° , L ( ” ) 

°H (2 )/ a L ( ”) 


a H (6) / T L ( " ) 

(eq. (37)) 

0.0 

0.000 

0.000 

0.000 

0.000 

0.000 

.2 

.057 

.220 

.215 

.217 

.212 

.5 

.118 

.266 

.271 

.271 

.268 

1.0 

.147 

.193 

.207 

.208 

.214 

2.0 

.118 

.117 

.123 

.127 

.132 

4.0 

.0693 

.0672 

.0673 

.0689 

.0713 

6.0 

.0478 

.0469 

.0466 

.0471 

.0483 


tThe first-order collision time r is related to t ct by 
the ratio ' r af T - 32(37r)“^y E , where y E is given in reference 6; 
in particular, y E (°°) yields 975 for the conversion of 

the independent variable in equation (37). 
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HIM 


TEMPERATURE EFFECTS 


Thermal Conduction 

Temperature gradients provide third and fourth checks of the collision model of 
equations (10) and (14) or (25) through comparisons with Spitzer and Harm's thermal con- 
ductivity and thermal diffusion. 

By a procedure similar to that used in appendix C, the Lor entz -like first-order 
Boltzmann equation for a plasma subject to a temperature gradient (but no pressure gra- 
dient) can be written 


Vf f 


(o) = f (o)^2 _ JjV. 


y 


R 


"2 Ce 


V In T P = -- 


13_ J-4/^ f (o) 


T o^-04 


<Pt 


(39) 


The solution of equation (39) is 


<Pf 


, \l/2 


which yields, respectively, the real and Lorentz (£ = 1) kinetic -energy fluxes 
q = \ n e m e ( c e 2 Ce) = -k 2 T e cr(e£) -2 ^ 2 + 5£ + 4)vT e 

and 

q L = -10k 2 T e e-2y E (oo)- 1 a vT e 

and, hence, the Spitzer-Harm thermal reduction factor 

■ 1 , 


5 t = (^l)' 1 ? = rE^ 00 ^ 10 ^ 2 ) (^ 2 + + 4 ; 


(40) 


(41) 


(42) 


(43) 


If the same value of | is substituted into equation (43) which was derived previously 
in equation (25) for the nonmagnetic zero -temperature- gradient problem, there results 


6 t (£ = 1.6674) = 0.3162 


(44) 


This time, however, the comparisons with Spitzer and Harm's exact value of 0.2252 and 
Kaneko's second Sonine approximation of 0.1906 are fairly poor. Even so, the value in 
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equation (44) is better than the corresponding Krook number of 0.0736 obtained from the 
substitution (see appendix B) of yg(l) and £, = 4 into equation (43). 

Thermal Diffusion 

Spitzer and Harm's thermal diffusion reduction factor also can be derived 

from equation (40) by first taking the velocity moment and then comparing real and 
Lorentz plasmas in the manner of equations (41) to (43). The pertinent expressions are 

T = ka(2e£) _1 (4 - £)VT e (45) 

y T = (T l ) _1 T = y E H(3^) _1 (4 - |) (46) 

and 


y T (! = 1.6674) = 0.2712 


(47) 


Despite the excellent agreement between the value in equation (47) and the exact 
value (ref. 6) of 0.2727, the preceding calculation is a valid test of equations (10) and (14) 
only in that it demonstrates the automatic incorporation by the present model of Onsager's 
reciprocal relations (ref. 9, pp. 704-717). This fact is illustrated quite directly by the 
following manipulation of equation (13) for no temperature gradient and of equations (41) 
and (45) for only a temperature gradient: 



q - h = 


ck2 Te (|2 + 5{ + 4, _ 5 


e 2 { 2 


v T t - g n e kT e v e 


ok 2 T (i; 2 + 5I+4) 5kT 

- L ! 


e2|2 


VT P + — — ^ j 
e 2e J 


2t 2 


ak^T 


( 9^2 + 16 ) 


4e £j 


„ „ V In Tp 

2^2 e 


(49) 
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The linear combination 


crkT e (4 - |) _ ok T p ; 
h = e 


2e^ 


? (9g 2 + 16) 


4e 2 | 2 


V In T e 


(50) 


is appropriate when the field and temperature gradients are present simultaneously. 

Likewise, the addition of equation (45) and Ohm's law for no temperature gradient 
gives the current density 


r - okT e (4 - « , 

j = ae + — gii — ln € 


(51) 


so that the matrix representation of equations (50) and (51) 



kT e (4 - 0 


\ 


2e^ 

k 2 T e 2 (9£ 2 + 16)1 


4e 2 £ 2 


/ 



(52) 


clearly reveals the symmetry between the off-diagonal elements. 

The fact that Onsager's reciprocal relations are guaranteed in the present proce- 
dure is, of course, another important advantage over several previous models. In par- 
ticular, the Krook model (appendix B), even though it technically satisfies the reciprocal 
requirements, does so in a trivial fashion because it predicts zero values for the off- 
diagonal matrix elements by forcing £ to equal 4. 


ALTERNATE DETERMINATION OF £ 


Equation (51) immediately suggests a more convenient experimental method for 
determining the value of £ in equation (10) than the method offered by equation (14). If 
a temperature gradient alone is initi ally applied to an enclosed plasma, a secondary gen- 
eralized electric field T will build up until the current produced by the temperature 
gradient is canceled (ref. 7, p. 144). The final stationary state so obtained corresponds 
to j = 0 in equation (51) and hence to the following expression for | in terms of the 
balancing e" and VT e : 


I = 



(53) 
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The values of £ obtained from equation (53) are, of course, numerically equivalent to 
those obtained from equation (14) because of the guaranteed satisfaction in the present 
technique of Onsager's reciprocal relations. 

As in equation (14), the conditions underlying equation (53) are assumed only for the 
purpose of establishing £, which remains fixed irrespective of subsequent conditions and 
forms of cp e if the plasma constituents do not change. It is difficult to imagine a sim- 
pler experimental procedure for completely defining the collision model than that of 
equation (53) and a corresponding one for the determination of t 0 in equation (10). The 
latter involves measurements of a field-induced electric current and the subsequent use 
of Ohm's law and equation (11). 


CONCLUDING REMARKS 

The collision model derived herein has been shown to predict quite adequately such 
diverse plasma properties as entropy, entropy and diffusion relaxation times, electro- 
magnetic tensor conductivities, and either thermal diffusion or off-diagonal kinetic-energy 
flux. The model thus provides a satisfactory scheme for utilizing macroscopic experi- 
mental data when microscopic parameters (e.g., electron — neutral-atom force laws) are 
scarcely known and for directly extending previous calculations to more complex condi- 
tions (e.g., to the fuli electromagnetic field). Only the simplest measurements, that is, 
electron diffusion and static field and temperature conditions, are required for the first 
application. In addition, the closed-form distribution functions are easier to interpret 
than the functions from polynomial expansions and reveal, for example, the natural 
appearance of entropy relaxation in Hall conductivities. 

Generalizations of this technique can apparently be made to higher perturbation 
orders. Certain modifications to the equations defining the collision model and to the 
kinetic equations are necessary, of course, but the manipulations should be much easier 
in a Lorentz-like framework than with the more exact Boltzmann collision integrals. 
Perhaps one or two second-order constraints will suffice to overcome Chapman-Enskog 
convergence difficulties just as the introduction of an effective collision parameter £ 
(empirical) overcame electron-electron difficulties in the present report. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., February 3, 1969, 

129-02-01-01-23. 
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APPENDIX A 


COLLISION INTEGRALS 


It is shown in this appendix that the last part of equation (6) follows directly from 
the second part when the electron — heavy-particle interaction potential can be written as 
the inverse |-power of the separation distance. 

The first step in such a development is the formulation of the reduced electron par- 
ticle velocity after a collision y T in terms of the reduced electron particle velocity 
before the collision y and the corresponding scattering or deflection angle x* This 
straightforward exercise in geometry yields 

y - = (1 - cos x)r - y sin x|^(sin cp cos e + cos 0 cos cp sin e) 

- j(cos cp cos e - cos 0 sin cp sin e) - k sin 0 sin e] (Al) 


where 0 and cp are the polar and azimuthal angles, respectively, of y in spherical 
coordinates and e is the azimuthal angle of y f with respect to the direction of y. 

Accordingly, equation (6) may be derived as follows with the aid of equation (5): 



The remaining integral involves, of course, the detailed collision dynamics, but that cal- 
culation has already been performed for the inverse power potential by Hirschfelder , 
Curtiss, and Bird (ref. 9, pp. 546-549). Their results can be written 


so that 



- cos x)t> db = 


/ X 1 / 2 

*U)/ m e \ -4/| 

27rnj\2kT e J 



= -K y 


1-4/ (o) 


<Pc 


(A3) 


(A4) 
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It must be mentioned, however, that the Coulomb potential (£ = 1) satisfies equa- 
tion (A3) only if a cut-off (usually the Debye length) replaces the upper limit of the inte- 
gral and if the y2 appearing in the argument of the logarithm in equation (8.3-8) of ref- 
erence 9 is replaced by its average value over all collisions. Both practices are rather 
standard in the current literature. 
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APPENDIX B 


THE KROOK COLLISION MODEL 


Holway (ref. 1) has used some concepts from information theory to derive the Krook 
collision model defined by 



m c 


- n, 


e l2T7kT, 


e h 


^3/2 r 

exp 


m. 



(Bl) 


More specifically, he obtained this result by means of a variational procedure based on 
the maximization of entropy subject to mass conservation and to experimental knowledge 
of the average velocity u e j and the average energy 3kT e j^/2 with which electrons 

emerge from encounters with their effective collision partners j. The pertinent rela- 
tions are equations (46) to (49) of reference 1, from which it is evident that T e j equals 
the electron temperature T e in the first Chapman-Enskog perturbation order. 

Other first-order expressions include 


f e - + <Pe) 

and the following linear expansion of equation (Bl): 



= -Kf e (0) ^ e - 2 ^ej 



(B2) 


(B3) 


where /3 e j and y are the reduced velocities 


,1/2 


m f 


e] " \2kTc 


^ei " 


*ej 


(B4) 


and 



(B5) 


The present objective is to determine 
equation 


<p e from the corresponding perturbation 


e 

m e 


E 


9f e <°) ( S£ e ^ 


= -Kf e (o) (<p e - 2^ e j * f) 


(B6) 
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for an applied electric field. If E is related to the reduced electron diffusion velocity 

a/2 


- / m. 


0 = 


>2kTV 


(B7) 


through Ohm's law 


— J m e -► m e _ 

E = ?= A^ , = '^ Ve '' 


m e / 2kT e 
eT o\ m e 


a/2 


(B8) 


there results 


1 _ 9f (o) 2f J°) „ _ 

4- p - n ‘ 


(o) 


'a dy 


t 0 • y = -Kl e ^'Up e - 2j e] ■ y) 


(B9) 


and 


’’e = 2 fejK + N • y 


(BIO) 


The definition 


'-®-M5V”Vedc e 


(°), 


(Bll) 


is used with equation (BIO) to yield 



and 


(B12) 


cp e = 2p-y (B13) 

and with equations (B3) and (B13) to yield 

= -ff e (0) <? e ( B14 ) 

a 

The physical interpretation of equations (B13) and (B14), and therefore the first- 
order Krook model, is clear: They correspond either to the first Sonine approximation to 

a real plasma or to the exact first-order solution for a Lorentz gas (no electron-electron 
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collisions) composed of Maxwellian molecules. Neither is adequate. Also in line with 
these conclusions is the observation from equation (10) that £ = 4 and = t ct (1) must 
be employed in the collision model of this paper if Krook model type of results are 
desired. 
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APPENDIX C 


CHAPMAN- ENSKOG SOLUTIONS 


The standard Chapman-Enskog first-order perturbation equation (ref. 2) for elec- 
trons can be written as follows for pressure gradients and an applied electromagnetic 
field E and B, no temperature gradients, infinitely massive heavy particles at rest 
relative to the laboratory, and the collision model of equation (7): 


'2kT 


1/2 


m. 


- y • Vf 


(o) e / m e 

e “ m e \2kT e/ 


vl/2 


9f 


(o) 


E 


by 


u>f e (o) (y x b) • = -Ky^^fg^Ve 


Since 


cp e = -jl^My + D 2 (y)y x B + D 3 (y)y • B&J • 

is a standard solution (ref. 3) of equation (Cl), there results 

. 1/2 


(Cl) 


(C2) 


— i 

m e \2kT ( 


^ ^ ^ y + w^D2y • BB - D 2 y + D jy X B^ = Ky^ - ^/^Diy + D 2 ~ X B + D 3 y • BB 


(C3) 


(o) 


upon taking the spatial and velocity derivatives of f e ' ' and defining the effective driving 
force as 


T = E + (en e )~^Vp 


(C4) 


It is readily confirmed by equating coefficients of y, y x B, and y • BB in equa- 
tion (C3) that the exact closed-form solution is 


cp e 


2 e / m e \ / 2 y/t)-lL + u£ 


m e K\ 2 kT e 


K 


^l + 4 y (8/?)-2 W i + Ti + « y (4/ ? )-lg 


K 


1 K 


B X £ 


(C5) 


where ej| and <T^ are the components of e" which are parallel and perpendicular, 
respectively, to B. The current density then becomes 
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j = -en e 




e2n e R 13 (~ ]13 r 

m e RR 04 \ 11 r 13 1 


wl 2 2 

+ KRi3 


B x e 


(C6) 


so that 


e2n e R 13 _ r 13 
m e CTR 04 t <j r 04 


(C7) 


Other definitions include the electrical conductivity a = e 2 n e T CT /m e and the integrals 




X-l 


1 + w2t ct 2 R 04 2 R 13 -^x 


(8/0-2)" >x2 ^ 


(C8) 


and = Iij( w = 0), 

Except for the determination of the following expressions thus represent the 
complete specification of the zero-temperature-gradient plasma: 


= - 


and 


2ctR 04 f m e 


en e R 13^kTp 


£ 


1/2 


(4 / D - l /, m 2 t o 2 r 04 2 ( 8/|)-2 

2 ^ 
r 13 Z 


-1 


' M 2t ? 2Rq 4 2 ( 8A) . 2 \ MT q?04 y (4/t)-lfe 

R 13 2 / " 1 R 13 


X 6 


(C9) 


-r U -1 _9 

j - a ( € || + r 13 *13 ^ + Ci,t ct r 04 R 13 i 22 b x 6 


(CIO) 


Also discussed in the main body of the report are the corresponding second Sonine 
approximation formulas derived from the matrix elements of Schweitzer and Mitchner 
(ref. 3). If the parameter a is introduced in these expressions for the purpose of 
describing the Lorentz (a = 0) and real ^a = 2^/2^ plasmas simultaneously, there results 


= 


4a 


m r 


en e (4a + 13)\2kT e 


<J( 3 y 2 + 2a - l)T 


+ 


-[ 3 (4a + 23)y 2 + 8a 2 + 22a - 43 B x 7> 
4a + 13 1_ 


(Cll) 
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T „ - . /l6a 2 + 104a + 259\£, v ~ 

j = a e + wr ff B x e 


\16a 2 + 104a + 169 


Only the linear terms have been retained in an coT a expansion. 
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THE ONSAGER MAXIMIZATION PRINCIPLE 

It is shown in this appendix that the perturbation function of equation (8) can be 
deduced from Onsager's principle on the maximization of the collisional production rate 
of entropy density (ref. 5). The mathematical statement of this principle is as follows: 



* kK6y y 1 " 4/4 f e (o) (Pe 2 e = 2kKyy 1-4/ ^f e (o V e 6<P e dc e = 0 (Dl) 

in which equation (7) is substituted for the collision integrals. 

If the diffusion velocity is taken to be the only nonequilibrium parameter on 
which cp e depends, the perturbation function can be varied arbitrarily in equation (Dl) 
subject to the constraint 

6j3 = &(y) = eng-^e^^l + <p e ) dc e = ng" 1 ^^ 0 ^ d^ e = 0 (D2) 

Accordingly, the application of the method of Lagrange multipliers to equations (Dl) and 
(D2) yields 

ty 1_4// ^f e ^6^ e ^C' e - ju • yy( 4 /^ -1 ^ dcfe = 0 (D3) 


and 


<Pe = • Y = 2R 04 R 13 “ V 4// ^ 1 0- Y (D4) 

Since equation (D4) is identical with equation (8), Onsager's maximization principle 
is embodied in the Boltzmann equation. It is further suggested that the present varia- 
tional procedure may be more efficient than the maximization of entropy used in informa- 
tion theory (ref. 1). At least equation (Dl) is based upon an objective physical theory, 
whereas the alternate technique refers to the uncertainty of the observer and provides 
only an upper limit to the true entropy unless the constraints form a complete set for 
describing the macroscopic state of the system. 
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THE MEAN- FREE -PATH MODEL 


The electron equation of motion (ref. 7, p. 28) in the mean-free-path statistical 
model can be written as follows for infinite ionic mass, an ionic charge number of unity, 
and zero-plasma-flow velocity: 


T= ct T - (en e ) -1 J X B - m e (e 2 n e ) 


-111 

0t 


(El) 


It is observed, however, from the macroscopic equations of continuity and energy transfer 
that the time derivatives of both particle density and temperature are second order in that 
they involve products of gradients and squares and divergences of the electron diffusion 
velocity; consequently, since the partial time derivative of equation (2) involves only 
0n e /0t and ST e /8t, there can be no time derivative on the left side of the first-order 
Chapman-Enskog perturbation equation (see appendix C). This argument means that the 
9 j j 9t appearing in reference 7 must be set to zero for comparisons with the present 
report. 

The cross product of equation (El) and B and the subsequent combination with 
equations (11) and (El) yield 

j x B = -ctBB x T + oB 2 (en e )”^B x f j x B 

= -en e u)T a B x 7 + en e co 2 T CT 2 (7^ - (en e )”* j* x b| (E2) 


the solution of which is 


Jx B = en e u>T CT |l + o) 2 t ct 2 ^ - Bxej 


A further combination of equation (E3) with equation (El) then gives 


j = o - 

and the tensor conductivities 


II + (l + co 2 r a 2 ) (7^ + ojt g B x ~ej 


(E3) 


(E4) 


a,. = a 


(E5a) 
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= CT 1 


+ W 2 r r 



ctjj = wr a a^l + a)^T a ^ 



Except for the substitution of r g for t q , equation (E4) is the same as equa- 
tion (35) to linear u>t q . Equation (E4) also gives the first Sonine approximation to 
when a assumes the corresponding value cr(l). 


(E5b) 

(E5c) 

T 
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